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Abstract 

Lorentz Invariance Violation (LIV) manifesting itself by energy dependent modification of stan- 
dard relativistic dispersion relation has recently attracted a considerable attention. Ellis et al. 
previously investigated the energy dependent time offsets in different energy bands on a sample of 
gamma ray bursts and, assuming standard cosmological model, they found a weak indication for 
redshift dependence of time delays suggestive of LIV. 

Going beyond the ACDM cosmology we extend this analysis considering also four alternative 
models of dark energy (quintessence with constant and variable equation of state, Chaplygin gas 
and brane- world cosmology). It turns out that the effect noticed by Ellis et al. is also present in 
those models and is the strongest for quintessence with variable equation of state. 
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I. INTRODUCTION 



Modern approaches to quantum gravity predict Lorentz Invariance Violation (LIV) as a 
consequence of the space-time foamy structure at small scales. Such violation can manifest 



0,6. 



itself by energy dependent modification of standard relativistic dispersion relation 

Several years ago Amelino-Camelia et al. [l| proposed to use astrophysical objects to look 
for energy dependent time of arrival delays. Specifically gamma ray bursts (GRBs) - very 



high energetic events visible : 
of constraining LIV theories 
LIV due to extra dimensions. 



rom cosmological distances - are the most promising sources 



6]. Recently Zhuk et al [7] have also used GRB bounds on 



The idea of searching for time of flight delays is tempered however by our ignorance 



concerning intrinsic delay (at the source frame) in different energy channels (see e.g. ^]). 
This clearly disfavors using energy dependent patterns in time-of-flights from single sources. 
As a possible way out of this trouble (and a first step in disentangling intrinsic time delays) 
Ellis et al. [4] proposed to work on a statistical ensamble of GRBs and formulated the 
problem in terms of linear regression where the intercept represents intrinsic time delay and 
the linear term represents LIV effect. They found a weak evidence for LIV (in [4j further 
corrected in [8]). On the other hand it was also shown {{j], that lack of detailed knowledge 
about cosmological model (in the context of accelerating expansion of the Universe) could 
be another source of systematic effects at high redshifts. Moreover, despite calling ACDM 
scenario "the concordance model" there exists voices - the most recent belonging to Linder 
[lo| - to consider several alternative scenarios such like the quintessence (scalar field with 
evolving equation of state) on equal footing in as many cosmological tests (CMBR, LSS, 
lensing etc.) as possible and only then judge which one is better supported by the data (see 



also 
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Hial) 



This suggests to ask how does such fitting of time delays vs. certain (defined below) 
function of cosmic expansion rate look in those alternative models. Or maybe the effect 
noticed by Ellis et al. could be attributed to incorrectly assuming the ACDM background 
whereas the real Universe's geometry is more accurately described by some other model (like 
the quintessence for example) taken from the inventory of accelerating cosmologies. Perhaps 
the redshift dependence of time delays vanishes in some other realistic (i.e. non-excluded by 
cosmological tests) models of presently accelerating Universe. This is the main idea behind 



the present paper. 



II. LIV INDUCED TIME DELAYS IN COSMOLOGICAL SOURCES 



Let us consider a phenomenological approach for LIV by assuming (after [2| for better 
comparision of results) the modified dispersion relation for photons in the form: 

E 



E l 
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where:e = ±1 is "sign parameter" 2j, £ n is a dimensionless parameter. One may assume (as a 
first guess) Eq G equal to the Planck energy £i = 1 and £ 2 = 10~ 7 [6j. The dispersion relation 
([[]) essentially corresponds to the power- law expansion (see [3j]) so for practical purposes 
(due to smallness of expansion parameter E/Eqg) only the lowest terms of the expansion 
are relevant. Because in some LIV theories the odd power terms might be forbidden [ijj] 
usually the cases of n = 1 and n = 2 are retained. For sake of being comparable with 
the results of Ellis et al. js] and because the LIV effects are small for sources at low and 
moderate redshifts used in this study we retain only the n — 1 term, and moreover we 
assume the sign parameter e = +1 in the formulae. 

Time of flight for the photon of energy E is equal to g, |8|, |9 ] 
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Consequently the time delay between a low energy and a high energy photon with the 
energy difference AE takes the following form 

AE r (1 + z')dz' 



At 
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where: H(z) = H h(z) is cosmological expansion rate (the Hubble function). 

The observational strategy emerging form ([3]) is very simple: monitor appropriate (i.e. 
emitting both low and high energy photons) cosmological source at different energy channels 



151 ] ) used a sample of 



and try to detect this time delay. Following this line Ellis et al. (j^, 4, 
gamma-ray bursts (GRBs) with known redshifts. However there remains an indispensable 
uncertainty: there is no reason for which low and high energy signal should be emitted 
simultaneously, and while detecting distinct signals (peaks in the light curve) at different 
energies we have no idea which one was sent first. This is known as so-called intrinsic time 
lags problem. 



In {4] it has been noticed that while this ambiguity clearly disfavors using energy depen- 
dent patterns in time-of-flights from single sources, one is still able to search for statistical 
correlations of spectral time lags with redshift in an ensamble of sources located at different 
redshifts. If one decomposes the observed time delay At D b s between different energy chan- 
nels: At obs = At LIV + At intrinsic then, using the notation of |4|, after taking into account 
cosmological time dilation factor 1 + z, one has At D b s = auvjj- + Z )K + b(l + z), where: 
K = and auv = HaE G ' ^ e °"S^ na ^ P a P er p| contained an error in the 

formula defining K(z) function. This was later corrected in j§| (correct expression was also 
given in |9j), but apart from fixing an obvious mistake the correction only changed numerical 
values of the results and was not able to erase the effect. 

Such parametrization allows to formulate the problem in terms of linear regression: 

Aids 



l + z 



a LIV K{z) + b (4) 



where the intercept informs about intrinsic time lags, and slope carries information about 
LIV effects. 



III. COSMOLOGICAL MODELS TESTED 



One of the most important issues in modern cosmology is the problem of "dark energy" , 
which appeared after the discovery of accelerated expansion of the Universe as inferred from 
the SNIa Hubble diagram If]]. Since then a lot of specific scenarios have been put forward as 
an explanation of this puzzling phenomenon. They fall into two broad categories: searching 
an explanation among hypothetical candidates for dark en ergy (cosmological constant A 
. quintessence - evolving scalar fields 17 1, Chaplygin gas 18 1) or modification of gravity 
theory (supergravity [3], brane world scenarios j^oj]). 

We will restrict our attention to flat models k = because the flat FRW geometry is 



strongly supported by cosmic microwave background radiation (CMBR) data 2l|,|22|. Fried- 
man - Robertson - Walker model with non-vanishing cosmological constant and pressure-less 
matter including the dark part of it responsible for flat rotation curves of galaxies (the co 
called ACDM model) is a standard reference point in modern cosmology. Sometimes it 
is referred to as a concordance model since it fits rather well to independent data (such 
like CMBR data, LSS considerations, supernovae data). The cosmological constant suffers 



from the fine tuning problem (being constant, why does it start dominating at the present 
epoch?) and from the enormous discrepancy between facts and expectations (assuming that 
A represents quantum-mechanical energy of the vacuum it should be 55 orders of magnitude 



larger than observed 



23|) 



Hence another popular explanation of the accelerating Universe is to assume the existence 
of a negative pressure component called dark energy. One can heuristically assume that 
this component is described by hydrodynamical energy-momentum tensor with (effective) 
cosmic equation of state: p = wp where — 1 < w < —1/3 [24|. In such case this component 
is called "quintessence". Confrontation with supernovae and CMBR data 25] led to the 



constraint w < —0.8. This was further improved by combined analysis of SNIa and large 
scale structure considerations (see e .g. 



26|) and from WMAP data on CMBR |2l|. The 



13| pinned down the equation of state parameter to 



ESSENCE supernova survey team 
the range w = —1.07 ± 0.09(stat) ± .12(systematics) . The most recent estimate of w 



— 0.969t o'^ll(stat)^Q^ll(systematics) comes from the Union08 compilation |27j. For the 



illustrative purposes we chose w = —0.87 as representing a quintessence model which is 
different enough from cosmological constant and still admissible by the data. 

If we think that the quintessence has its origins in the evolving scalar field, it would 
be natural to expect that w coefficient should vary in time, i.e. w = w(z). An arbitrary 
function w(z) can be Taylor expanded. Then, bearing in mind that both SNIa surveys 
or strong gravitational lensing systems are able to probe the range of small and moderate 
redshifts it is sufficient to explore first the linear order of this expansion. Such possibility, i.e. 
w(z) = w + W\z has been considered in the literature (e.g. 1 28]). Fits to supernovae data 
performed in the literature suggest wo = —1.5 and w\ = 2.1 [29] (which is consistent with fits 
given in [l2|). Therefore we adapted these values as representative for this parametrization 
of the equation of state. 

In the class of generalized Chaplygin gas models matter content of the Universe consists 
of pressure-less gas with energy density p m representing baryonic plus cold dark matter 
(CDM) and of the generalized Chaplygin gas with the equation of state pch = — with 
< a < 1, representing dark energy responsible for acceleration of the Universe. Using the 
angular size statistics for extragalactic sources combined with SNIa data it was found in 30] 
that in the the fl m = 0.3 and Qch — 0.7 scenario best fitted values of model parameters are 
A Q = 0.83 and a = 1. respectively. Generalized Chaplygin gas models have been intensively 



studied in the literature 
data (e.g. 



31] and in particular they have been tested against supernovae 



32| and references therein). Conclusions from these fits are in agreement with 



the above mentioned values of parameters so we used them as representative of Chaplygin 
Gas models. 



Brane-world scenarios assume that our four-dimensional spacetime is embedded into 
5-dimensional space and gravity in 5-dimensions is governed by the usual 5-dimensional 
Einstein-Hilbert action. The bulk metric induces a 4-dimensional metric on the brane. The 
brane induced gravity models 20[ have a 4-dimensional Einstein-Hilbert action on the brane 
calculated with induced metric. According to this picture, our 4-dimensional Universe is a 
surface (a brane) embedded into a higher dimensional bulk space-time in which gravity prop- 
agates. As a consequence there exists a certain cross-over scale r c above which an observer 
will detect higher dimensional effects. Cosmological models in brane-world scenarios have 
been widely discussed in the literature [33JI. It has been shown in 33] that flat brane-world 
Universe with Vt m = 0.3 and r c = 1.4 Hq 1 is consistent with current SNIa and CMBR 
data. Note that in flat (i.e. k — 0.) brane-world Universe the following relation is valid: 
^r c — z(l — ^m) 2 - Futher research performed in |34j based on SNLS combined with SDSS 



35] using 



disfavored flat brane-world models. More recent analysis by the same authors 
also ESSENCE supernovae sample and CMB acoustic peaks lead to the conclusion that 
flat brane-world scenario is only slightly disfavored, although inclusion of baryon acoustic 
oscillation peak would ruled it out. Despite this interesting debate we use flat brane-world 
scenario with Q m = 0.3 for illustration. 



Because in what follows we will use model selection approach, it would be tempting 
- indeed even more appropriate from the methodological point of view — to take model 
parameters as free and estimate their values from the best fits. However, with a small 
sample and considerable parameter space, it is almost sure that the values (e.g. Q m , w, 
Wo, W\, etc.) inferred would be inconsistent with general knowledge already acquired with 
alternative techniques on much better data and much larger samples (like SNe la, CMBR, 
LSS). The major strength of modern cosmology is in consistency, therefore we take fixed 
values (best fits to SN la data) of cosmological parameters specifying the models considered. 

6 



IV. SEARCHING FOR LIV EFFECTS IN DIFFERENT COSMOLOGICAL MOD- 
ELS 



The sample we used consists of 35 GRBs with known redshifts for which time lags between 
different energy channels have been assessed from the light curves by Ellis et al. ^]. They 
are summarized in Table 1 of [4]. The data are based on the results of BATSE (9), HETE 
(15) and Swift (11) experiments. The numbers in brackets indicate the number of bursts 
detected in a specific experiment. Technical details can be found in ^|. We took these data 
for the sake of comparability. 

The dependence on cosmological model resides in the K(z) term which was calculated 
for all models considered. Their expansion rates expressed through cosmological model 
parameters are shown in Table I. The parameter values we used were described in previous 
section. The results of linear regression are summarized in Table II and displayed in Figure 
1. One can see that the LIV effect (i.e. non-zero slope) is visible in all models. For most 
of them it is significant at 2a level and for the quintessence model with varying equation of 
state at 3a level. One should note, however, that our regression is a straightforward one not 
supported with sophisticated tricks original authors made after such straightforward fit in 
order to get their final estimate. Therefore, we do not attempt to convert slope factors into 
LIV energy scale estimate. 

One can see that in all classes of alternative cosmological models the effect is similar, 
hence there is no indication that Ellis et al. result might be an artifact of assuming ACDM. 
Now, one can ask which cosmological model "is the best" in the sense of revealing a sought- 
for LIV signature i.e. the linear relation (j3J). 

In order to compare different models - the problem we encounter here - one can use 
Akaike information-theoretical model selection criterion 137 
become a standard diagnostic tool of regression models [38 



applied by Liddle 



39j and then was used in e.g. 
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40 



4l|- 



. In particular this criterion has 
. In cosmology it has first been 



Akaike criterion is based on Kullback-Leibler information /(/, g) between two distribu- 
tions f(x) and g(x). The intuitive meaning of I(f,g) (also called K-L divergence) is the 
information lost when g is used to approximate /. It is convenient to think that f(x) denotes 
the true mechanism behind the data and g(x\9) its approximating model (parametrized by 
6). Of course, K-L divergence cannot be assessed without prior knowledge of the true model 



f(x) as well as parameters 9 of the approximating model g(x\6). However, given f(x) and 
g(x\8) there exists the "best" value of 6 for which Kullback-Leibler divergence is minimized. 
The maximum likelihood estimator 8 of 6 parameter is exactly this K-L "best" one. 

The core result of Akaike was in showing that an approximately unbiased estimator of 
K-L divergence is ln(C(9\data)) — K where £ is the likelihood function (more precisely its 
numerical maximum value - taken at 6) and K is the number of estimable parameters (9) 
in approximating model g(x\8). For historical reasons Akaike formulated this result in the 
following form: 

AIC = -2ln(C(9\data)) + 2K (5) 

which became known as Akaike information criterion. Heuristically one may think of it as 
of an estimator of K-L divergence between the model at hand g(x\8) and an unknown true 
model f(x) which generated the data. The first term measures goodness of model fit (or 
more precisely the lack thereof) and the second one (competing with the first) measures 
model complexity (number of free parameters). 

The AIC value for a single model is meaningless (simply because the true model f(x) 
is unknown). What is useful, instead are the differences Aj := AICi — AIC min calculated 
over the whole set of alternative candidate models % = 1, N where by AIC m i n we denoted 
min{AICi] i = 1, N}. Comparing several models, the one which minimizes AIC could be 
considered the best. Then the relative strength of evidence for each model can be calculated 
as the likelihood of the model given the data C(gi\data) cx exp(-^Ai). Relative likelihoods of 
the models C(gi\data) normalized to unity are called Akaike weights u>j. In Bayesian language 
Akaike weight corresponds to the posterior probability of a model (under assumption of 
equal prior probabilities). The (relative) evidence for the models can also be judged by the 
evidence ratios of model pairs — = ^ttt^- 

1 Wj L(gj\data) 



A very similar criterion was derived by Schwarz [42 
as the so called Bayesian information criterion (BIC) 



in a Bayesian context. It is known 



BIC = -2ln{C0\data)) + Kln{n) (6) 

where n is sample size and as previously K denotes number of parameters. BIC is not an 
estimator of K-L divergence - its derivation stems from estimating the marginal likelihood 
of the data (marginalized over parameters). BIC does not take the full advantage offered 

8 



by Bayesian techniques such as described e.g in Trotta and Kunz [43j (see also references 
therein) . 

At last, it should be noted that AIC may perform poorly if the sample size is not large 
enough (n < 40 as a rule of thumb), which is the case of our study. In order to cope 
with such situations Sugiura derived a variant of Akaike criterion called c-AIC, which 
is related with AIC in the following manner: AIC C = AIC + 2K( " I t~ l ) . 

In our case i.e. simple linear (univariate) regression the distinction between AIC, c-AIC 
and BIC is purely formal. Technically the corresponding numerical values are shifted by a 
constant, which does not matter since there are the differences that count. If one used the 
freedom of cosmological model parameters entering K(z) while performing the fit, clearly 
the c-AIC should be used with the sample size of n = 35. In our specific case of LIV 
testing this is not an option. Hence the reported AIC analysis summarized in Table III is 
representative for all above mentioned model selection criteria (and due to simplicity also 
to the standard goodness of fit techniques). The likelihoods needed to calculate AIC were 
taken here as exponents of minus chi-square for the best fitted model. One can see that 
the quintessence model with varying equation of state gives the best fit in time delay vs. 
K{z) regression. However the odds against the other models are not big enough to call it 
preferred. 



V. DISCUSSION AND CONCLUSIONS 

In this paper we performed fits of time delays (between different energy channels in a 
sample of GRBs) versus the K(z) function (dependent on the cosmic expansion rate) in five 
cosmological models representative to different dark energy scenarios. The technique we 



used was taken over an original one of Ellis et al. [4] and so was also the sample adopted. 
This was motivated by suggestion formulated in |4j that the LIV effect shows up in the 
studied sample of gamma ray bursts. These authors have used the "concordance" model 
(i.e. flat ACDM with = 0.3) as representing the cosmological background. 

Despite the fact that such ACDM model is commonly used as a reference point in cosmo- 
logical studies, there is not so obvious that this is the true solution of accelerating Universe 



enigma. See e.g. the most recent discussion of this issue by Linder 10j (and references 



therein). There exist several classes of alternative dark energy scenarios, from which we 



Table I. Expansion rates H{z) in the models tested. The quantities fij represent fractions of critical 
density currently contained in energy densities of respective components (like clumped pressure-less 
matter, A, quintessence, Chaplygin gas or brane effects). 



Model 



ACDM 
Quintessence 
Var Quintessence 
Chaplygin Gas 

Braneworld 



Cosmological expansion rate H{z) (the Hubble function) 



h\z) = hI [n m (i + z ) 3 + n A } 
h 2 (z) = Hi [n m (i + zf + n Q (i + z) 3 d+-)] 

H 2 (z) = Hi [Q m (1 + zf + £Iq (1 + z )3(i-H«o-i»i) exp(3^iz)] 



H{zf 



Hi 



H(z) 



n m (i + zf + n Ch (Ao + (i - Aq)(i + ^) 3(1+Q) ) 1+a 



Hi 



have taken four most commonly considered. 

Initially we suspected that the effect noticed in [4] could be an artifact of incorrectly 
assuming ACDM in their fits whereas another model better represents the accelerating Uni- 
verse. Contrary to our expectations the effect does not get smaller in alternative models. In 
fact it is the highest in the model providing the best fit of time delay vs. K(z) regression. 

Curiously the quintessence model with varying equation of state is the one which gives 
the best fit. It is most probably a mere coincidence but it is this cosmological model which is 
by many considered to correctly reflect the phenomenology of the dark energy. Dark energy 
is the solution of presently accelerating expansion of the Universe, and, as brilliantly noticed 
in [lO I , the only accelerating expansion phase we know (i.e. the inflation) has ended hence 
clearly had time variation and dynamics. 

Our approach was a very simple one, aimed mainly to raise the issue. However the results 
suggest that it would be interesting to go beyond taking compilation of and re-analyze 
for the LIV effects, assuming alternative cosmological models of an accelerating Universe, on 
the raw data (GRBs light curves in different energy channels) and with all the sophistication 
used in the original study. 
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Table II. Regression coefficients (with la ranges) for the time delay vs. K(z) technique in the 
cosmological models tested. 



Cosmological model 


Regression coefficient clliv 


Intercept b 


ACDM 
Quintessence 
Var Quintessence 
Chaplygin Gas 
Braneworld 


a LIV = -0.0794 ± 0.0447 
a LIV = -0.0806 ± 0.0460 
a LIV = -0.1510 ±0.0683 
a LIV = -0.1201 ±0.0618 
a LIV = -0.0866 ± 0.0493 


b = 0.0494 ± 0.0288 
b = 0.0489 ± 0.0288 
b = 0.0735 ± 0.0340 
b = 0.0627 ± 0.0330 
b = 0.0501 ± 0.0294 



Figure 1. Results of the linear regression for time delay vs K(z) technique in the cosmological 
models tested (the case of Quintessence was too similar to the ACDM to deserve displaying) . 
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